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SUMMARY 


Because of the increasingly complex environments in which the flight crews of 
commercial aviation aircraft must operate, a research effort is currently 
underway at NASA Langley Research Center to investigate the potential benefits 
of intelligent cockpit aids, and to establish guidelines for the application 
of artificial intelligence techniques to advanced flight management 
concepts. The segment of this research area that concentrates on automated 
fault monitoring and diagnosis requires that a reference frame exist, against 
which the current state of the aircraft may be compared to determine the 
existence of a fault. This paper describes a computer program which generates 
the portion of that reference frame that specifies the horizontal flight 
route. 


INTRODUCTION 

In the commercial aviation industry, safety and efficiency are of utmost 
importance. Flight crews must operate in increasingly complex environments, 
and automation in the form of intelligent cockpit aids may assist them in 
improving safety and efficiency. Modern microprocessor and display technology 
have made it feasible to automate many flight deck functions, and the advent 
of expert systems technology has proven that artificial intelligence (AI) has 
a great deal to offer to enhance and aid the man in the system. In 
particular, automation of fault monitoring and diagnosis is desirable to 
ensure safe and efficient operations, since response to faults on board an 
aircraft is often time critical. It is believed that the use of artificial 
intelligence concepts can facilitate this automation. 

An effort is currently underway at NASA Langley Research Center to investigate 
the potential benefits of intelligent cockpit aids, and to establish 
guidelines for the application of AI techniques to advanced flight management 
concepts. One application area within this research effort is onboard fault 
monitoring and diagnosis. The segment of this area that concentrates on 
automated fault monitoring requires that the actual aircraft state be compared 
to some desired reference aircraft state in order to detect the existence of a 
fault. 

The aircraft state has been divided into portions to facilitate the generation 
of the desired reference aircraft state. A computer program has been 
developed to generate the portion of this reference state that specifies a 
horizontal route minimizing ground distance between the origin and destination 
airports using the established airway system designated by the Federal 
Aviation Administration and in use by Air Traffic Control. This report 
documents this computer program and describes its use and capabilities. 


BACKGROUND 

The initial concepts for the fault monitoring and diagnosis segment of the 
intelligent flight management research were developed at the University of 


Illinois at Urbana-Champaign under a grant from NASA Langley Research Center 
(ref. 1). The major thrust of the work was to explore the existing artificial 
intelligence techniques applicable to fault monitoring and diagnosis within 
the flight domain and, if necessary, develop new techniques more suited to 
this application. Out of this research came the basic concepts for a multi- 
level architecture within which the fault monitoring segment would operate. 

The description and precise function of each of the levels within this 
architecture have evolved as new ideas were incorporated. 

The basic architecture of the fault monitoring segment is organized as a 
hierarchy. The hierarchy consists of four distinct but interconnected levels 
as shown in Figure 1: the Route Level, the Trajectory Level, the Flight 

Dynamics Level, and the Subsystems Level. In general, each level is 
subservient to the one above it, supplying the means of achieving the goals 
which are passed to it from above. This holds true except in the case of the 
Subsystems Level; the relationship here could be more accurately described as 
a precondition-satisfaction check. An expert system is included within this 
framework to handle all diagnostic duties. The system architecture, the 
monitoring functions, and the expert system are written in the LISP language 
dialect known as MACLISP, and the Route and Trajectory Level Generators are 
written in FORTRAN 77. The following is a general description of each level 
within the fault monitoring and diagnosis segment of the intelligent flight 
management research area. 


The Route Level 


The Route Level is the highest level in the architecture, passing along its 
goals and restrictions to the Trajectory Level. These primary flight goals 
include the generated route and any time, weather, traffic, or other Air 
Traffic Control (ATC) restrictions. The Route Level consists of two main 
sections, the Route Generator and the Route Monitor. The Route Generator 
generates the shortest route from the city of departure to the city of 
arrival, taking into consideration constraints such as weather conditions and 
ATC restrictions. The Route Monitor continuously monitors the progress of the 
flight to determine if the generated route is being followed. 


The Trajectory Level 

This level also consists of two sections, the Trajectory Generator and the 
Trajectory Monitor. The Trajectory Generator utilizes the primary flight 
goals sent down from the Route Level Generator. Once these are received, the 
Trajectory Generator produces an optimal trajectory that minimizes either fuel 
or direct operating costs. The generated trajectory is passed down to the 
Flight Dynamics Level. The Trajectory Monitor monitors the progress of the 
flight from the vertical- and time-based viewpoint, deciding whether or not 
the optimal trajectory is being properly followed. 
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The Flight Dynamics Level 


This level uses an internal model of the aircraft to generate all necessary 
control inputs such as wheel and column positions, throttle positions, and 
flap settings that are required to achieve the desired trajectory. This 
sequence of control inputs is then transferred down to the Subsystems Level. 


The Subsystems Level 

The Subsystems Level checks the control input sequence and amplitude to assure 
that the inputs do not attempt to extend the aircraft beyond its normal 
operational and performance limits. This level contains the knowledge about 
the aircraft subsystems, and checks that each subsystem is working properly. 

When the checks at the Subsystem Level have been completed, the Flight 
Dynamics Level passes its sequence of control inputs back up into the 
Trajectory Level, which then passes the generated optimum trajectory back up 
to the Route Level. Here it is stored for referencing by the route monitoring 
functions. Once the flight has been initiated, constant surveillance of all 
critical parameters is maintained, and if any condition is observed that has 
caused or will cause deviation from the planned route, the expert system is 
activated. The expert system itself currently uses a rule-based approach to 
problem solving as described in ref. 2, but this will later be expanded to 
more closely combine the functions of error detection and fault diagnosis into 
a more unified effort. 


DESCRIPTION OF THE ROUTE GENERATOR 
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current information on weather, traffic, and Air Traffic Control (ATC) 
instructions can be used. If any of these initial conditions or instructions 
are changed during the flight so as to render the previously generated route 
unacceptable, the Route Generator can be executed again using the next 
waypoint of the path as the point of departure, and the destination airport, 
whether different or the same as originally planned, as the point of 
arrival. This inflight re-routing capability is the principal advantage of a 
program which actively generates the route rather than retrieving a rigid, 
preplanned course from a storage device. The Route Generator itself is 
comprised of two major sections, the Selector and the Compiler. 


The Route Selector 


The Route Selector is the section which chooses the shortest allowable path 
from origin to destination using the established High Altitude Jet Route 
System. For the purposes of analysis, the Jet Route System is represented 
within this program in graph form, the nodes representing geographical 
reference points, or waypoints, and the arcs between these nodes denoting the 
airway segments. Each arc must be assigned a weighting factor. For the case 
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of locating the shortest path between any two nodes in the graph, this 
weighting factor is represented by the mileage between waypoints. The matrix 
containing the distance information is arranged so that arcs are identified by 
their location within the two-dimensional graph. The index for each dimension 
is a pointer to the name of the corresponding node (waypoint) in the one- 
dimensional matrix containing the node identifiers. Thus the graph, when 
printed out, takes the form of a mileage chart, with only those mileages shown 
that represent existing segments. A sample portion of the graph is shown in 
figure 2. The subroutine which performs the actual selection of the sequence 
of arcs from origin to destination nodes uses Dijkstra's algorithm (ref. 3). 
This algorithm locates the path of least resistance through an array of 
weighting factors, and is ideally suited for application to this particular 
class of problem. 


The Route Compiler 


The Route Compiler assembles the sequence of track angles (no-wind true 
headings) and the corresponding distances required to follow the chosen path 
through the established system of jet airways, as well as the magnetic 
variation from true north associated with each waypoint in the path. This 
sequence contains the route information needed by the Trajectory Level for 
generation of the optimum trajectory, since it describes the actual track of 
the aircraft across the ground. 


CONSTRUCTION OF DATA BASES 

There are four data bases accessible by the Route Generator. The first data 
base is contained in the matrix called NAME, which holds the names of all the 
waypoints stored in alphabetical order. The position of each waypoint name 
within this array is very important, because the remaining data bases use this 
position as a pointer for identifying segments between waypoints within the 
connectivity matrix. The second data base, called the primary data base, 
contains all the information necessary to permit the Selector to choose the 
shortest route through the airway system from origin to destination. This 
includes the endpoints of each segment, the segment length, and the directions 
of travel permissible on each segment. The third data base contains the track 
information necessary for the Compiler to assemble the entire sequence of 
track angles and mileages for passage down to the Trajectory Level. The 
fourth data base contains the magnetic variation from true north associated 
with each waypoint. 

Each entry in the primary data base lists general information on one 
particular airway segment. The entry is constructed entirely of integer 
numbers, and has the following format: 

AAABBB CCC ODD E 
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where 


AAA is a pointer to the location of the name of the first endpoint in 
the NAME array, 

BBB is a pointer to the location of the name of the second endpoint in 
the NAME array, 

CCC is the Jet Route designation number, 

DDD is the length, in nautical miles, of that entire segment, and 

E is a direction constraint flag indicating the allowable directions 

of travel on that segment: 

0 - open (travel in both directions permissible) 

1 - travel inbound to first endpoint only 

2 - travel inbound to second endpoint only 

For example, notice in figure 3 that Fort Lauderdale DME is named FLL. In 

this implementation FLL is located in the NAME array in position 21. Orlando 

VORTAC is named ORL, and is located in the NAME array in position 48. The 
segment connecting these two waypoints is designated J20, and has an overall 
length of 161 nautical miles. Therefore, this segment is stored in the 
primary data base as: 

021048 020 161 0 
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segment is normally open to travel in both directions. 


Each entry in the third data base also lists information on only one airway 
segment, and is constructed of only integer numbers, but here it has the 
following format: 


AAABBB FF GGGHHH IIIJJJ KKKLLL ... etc. 

where 

AAA is a pointer to the location of the name of the first endpoint in 
the NAME array, 

BBB is a pointer to the location of the name of the second endpoint in 
the NAME array, 

FF is the number of different portions comprising the airway segment, 
GGG is the track angle of the first portion. 
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HHH is the distance, in nautical miles, that the first track angle is 
to be followed 

III is the track angle of the second portion, and 

JJJ is the distance, in nautical miles, that the second track angle is 
to be followed . . . etc. 


Also notice in figure 3 the path J20 takes between FLL and ORL. J20 first 
tracks 339 degrees for 92 nautical miles, and then tracks 334 degrees for 69 
more nautical miles. This information for J20 is represented in the secondary 
data base as: 

021048 02 339092 334069 

Since the exact location of all waypoints and turning points must be known in 
order to calculate the track angles, the information for this data base must 
be gathered from other sources in addition to the Enroute High Altitude Charts 
which display the Jet Route System. The data base containing the magnetic 
variation from true north at each waypoint is stored in an array. Each 
element in the array has associated with it a VORTAC name and a character 
which is "W" to indicate a westerly magnetic variation and "E" for an easterly 
variation. 

The designated airways of the Jet Route System are not the only available 
paths through the sky; however, they are the most desirable ones from the 
viewpoint of ATC. Special clearances are routinely granted for flights to 
cities not served by the System, and for flights between cities where use of 
the System would obviously present a hindrance to both the airline flight crew 
and the ATC staff. In commercial applications of this flight planning 
program, therefore, the data bases will need to be tailored to the existing 
structure of allowable routing for each individual airline. 


RESULTS AND DISCUSSION 

Several sample runs were made to demonstrate the capability of the program to 
produce the shortest route between waypoints. Various features of the route 
generator program are illustrated, including the capability of reversing 
direction and of specifying that certain flight segments are undesirable. 

The first example illustrates a case which was run for the trip between Miami, 
Florida and Atlanta, Georgia. Figure 4 shows the interactive dialogue of the 
user with the route generator program and the output of the program. The user 
specifies the departure and arrival points and any undesirable segments. The 
program produces a list of the airways segments and corresponding mileages 
that make up the generated route. The waypoints listed include Miami VORTAC 
(MIA), LaBelle VORTAC (LBV), Lakeland VORTAC (LAL), and Atlanta VORTAC 
(ATL). The Jet Route segments listed include J43 between MIA and LBV, J73 
between LBV and LAL, and J89 between LAL and ATL. 
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The second example is the reversal of the first example trip. The departure 
point is Atlanta and the destination is Miami. Figure 5 shows the results of 
this example, which shows the reversing function of the subroutine GNRATR. 
GNRATR adds 180 degrees to the track angles and reverse their order whenever 
travel occurs in the direction opposite from that in which the track and 
mileage information is stored. 

Two sample runs that illustrate the ability of the program to plan around 
route segments that have been flagged as unacceptable are illustrated in 
figures 6 and 7. These runs demonstrate the fact that the program does always 
choose the shortest path between waypoints. Both cases were run for the trip 

between Hobby, Texas and Oklahoma City, Oklahoma. The waypoints listed 

include Hobby (HUB), Humble (IAH), Navasota (TNV), Dallas - Ft. Worth (DFW), 
and Oklahoma City (OKC). The Jet Route segments listed include J177 between 
HUB and IAH, J87 between IAH and TNV, J33 between IAH and DFW, 087 between TNV 

and DFW, and J21 between DFW and OKC. 


CONCLUDING REMARKS 

The computer program for generating a horizontal route that minimizes ground 
distance has been identified as part of the fault monitoring function of an 
onboard fault monitoring and diagnosis system. Using a small subset of the 
geographical reference points in the High Altitude Jet Route System, the 
program has been demonstrated to generate the shortest route between specified 
source and destination airports. The capability of the route generator 
program to route around path segments designated as undesirable was also 
demonstrated. 
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APPENDIX 


Description, Recommendations, and Source Listing of Program 
ROUTEGEN and Subroutines LOADER, NOGOOD, SLECTR, GNRATR , and MAG VAR 


The computer program described in this paper was written in FORTRAN 77 on a 
Prime 750 mini -computer, running under the Primos* operating system. A 
Tektronix 4014-1 was used for the remote terminal. The primary data base is 
stored on disk as the file named SEGDAT, and the secondary data base is stored 
as the file named TRACKDAT. 


Program ROUTEGEN 


Description 

This is the driver program for the subroutines which perform the actual 
functions of selection and generation of the route. The origin and 
destination are read in from this level, and the selected route is output from 
this level. 

Recommendations 


1. The program currently communicates interactively with the user via a 
CRT terminal. If desired, points of departure and arrival could be 
pre-programmed or read in from another source, and the output could 
be printed out on a NAV display or stored for future reference by the 
pilot. 

2. The call to the subroutine LOADER could be eliminated by storing all 
of the primary data base information in the appropriate arrays 
through the use of DATA statements. The execution time would 
definitely decrease, but the program size would increase. Depending 
on the overall size of the graph network, this may or may not be a 
desirable alternative. 

Source Listing 

The listing of program ROUTEGEN is as follows: 
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READ( IN, 20 )NAME1, NAMES 
FORMAT (A3, IX, A3 ) 
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END 


APPENDIX 


Subroutine LOADER 


Description 

This subroutine stores the information from the primary data base in the 
appropriate arrays. To minimize storage requirments, this information is 
given in the primary data base for travel in one direction only. That is, all 
the segments are stored heading generally south-to-north. As the information 
is being read into the different matrices, it is stored in both the south-to- 
north and the north-to-south positions. 

Recommendations 


1. The LENGTH array needs to be zeroed before the mileage information 
can be read in, and it has been done here using a DATA statement. 
This has the effect of increasing the storage requirements 
dramatically. A slower but more space efficient method would be to 
EQUIVALENCE the LENGTH array to a one-dimensional array, and then 
zero it using a single DO loop. If a machine-dependent function is 
available for this purpose it should be used in any implementation 
which still makes use of the LOADER subroutine. 

Source Listing 

The listing of subroutine LOADER is as follows: 
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-THIS SUBROUTINE IS USED TO LOAD THE SPARSE LENGTH MATRIX USED 
-TO FIND THE SHORTEST ROUTE BETWEEN SOURCE AND SINK 
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APPENDIX 


Subroutine NOGOOD 


Description 

This subroutine allows certain airway segments to be selected and labelled as 
unacceptable, for reasons of traffic, weather, or ATC constraints. The 
subroutine queries the user for the endpoints of each bad segment, and then 
inserts a zero for the mileage of that segment in the LENGTH matrix. This 
effectively removes the segment from consideration within the SLECTR 
subroutine. 

Recommendations 


1. All the information needed for labelling unacceptable segments could, 
if desired, be read in from a source other than interactive 
communication with a user. A realistic implementation of this 
program might have the traffic, weather, and ATC constraints stored 
in a dynamic data base, possibly being regularly updated by a data 
base monitor/modifier program. The Route Level Monitor would then 
need to periodically check the status of all segments in the planned 
route to determine if modification of the chosen flight plan is 
necessary. 

Source Listing 

The listing of subroutine NOGOOD is as follows: 
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APPENDIX 


Subroutine SLECTR 


Description 

This subroutine uses Dijkstra's algorithm (ref. 3) to find the path of least 
resistance (smallest total mileage) through the LENGTH array. Minor 
modification was necessary to incorporate the direction constraint flag 
capability, but the algorithm itself remains basically unaltered. 

Source Listing 

The listing of subroutine SLECTR is as follows: 
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SUBROUTINE SLECTR 

-THIS SUBROUTINE PICKS THE SHORTEST PATH THROUGH THE LENGTH MATRIX 
-AND STORES IT IN THE 'PATH' ARRAY. 
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-FIND THE TENTATIUELV LABELLED NODE UITH THE SMALLEST LABEL 
PIIN» INFINITY 


IF << STATE ( I, LABL KEG. TENT). AND. (STATE (I, LENGTH >.LT. MIN) >THEN 
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APPENDIX 


Subroutine GNRATR 


Description 


This subroutine compiles the track angle and mileage sequence from the third 
data base. The current technique used for this purpose employs a simple, but 
slow, process of searching the entire data base for each list of segment 
information. As each segment list is located, it is printed out to the CRT 
screen according to the direction of travel. If travel occurs in the 
direction opposite that in which the segment information is stored, 180 
degrees are added to the track angles, which are output in reverse order with 
their mileage. 

Recommendations 


1. A much faster search method would include storing the airway segment 
data in the same sequence within both data bases. An array of 
pointers to the location of each chosen segment could be filled at 
the same time as the array of pointers to the waypoint names, known 
as the PATH array. The track and mileage information could be 
compiled much faster with the few number of disk file operations 
required. Once assembled, the subroutine GNRATR currently merely 
prints the track and mileage information to the CRT terminal screen 
and to a file for later use by the Trajectory Level. When combined 
with other computer programs that perform the functions of the other 
levels in the fault monitoring and diagnosis segment's basic 
architecture, this information will be used by the Trajectory Level. 

Source Listing 

The listing of subroutine GNRATR is as follows: 
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SUBROUTINE GNRATR 

-THIS SUBROUTINE RETRIEVES FROH THE DATA BASE THE ENTIRE SEQUENCE 
-OF TRACK ANGLES (NO WIND TRUE HEADINGS) AND CORRESPONDING 
-DISTANCES THAT IS NEEDED TO TRAVEL ALONG THE PATH PREVIOUSLY 
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APPENDIX 


Subroutine MAGVAR 


Description 

This subroutine retrieves the magnetic variation of each point along the path 
from true north. The data base is searched using a simple sequential search 
until the desired waypoint names are found. The magnetic variation for each 
waypoint in the chosen path is stored in an array called VAR_MAG. 

Source Listing 

The listing of MAGVAR is as follows: 
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RETURN 

100 URITE ( 6, 101 ) 

101 FORMAT ( IX, '+•*■+ EOF FOUND IN UARIATION DATA FILE + + + 
RETURN 

END 
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Figure 1.- Multi-level architecture of the fault monitor. 
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Figure 2.- Partial printout of the contents of the graph matrix, 

illustrating the connectivity concept between waypoints. 
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Figure 4.- Example dialogue and output for the trip from Miami, Florida 
to Atlanta, Georgia. 
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Figure 5.- Example dialogue and output for the trip from Atlanta, Georgia 
to Miami, Florida. 
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Figure 6.- Example dialogue and output for the trip from Hobby, Texas 
to Oklahoma City, Oklahoma. 
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Figure 7.- Example dialogue and output for the trip from Hobby, 
Texas to Oklahoma City Oklahoma with undesirable 
segments specified. 
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